Superfluid, solid, and supersolid phases of dipolar bosons in a quasi-one-dimensional 

optical lattice 
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We discuss a model of dipolar bosons trapped in a weakly coupled planar array of one-dimensional 
tubes. We consider the situation where the dipolar moments are aligned by an external field, 
and find a rich phase diagram as a function of the angle of this field exhibiting quantum phase 
transitions between solid, superfluid and supersolid phases. In the low energy limit, the model 
turns out to be identical to one describing quasi-one-dimensional superconductivity in condensed 
matter systems. This opens the possibility of using bosons as a quantum analogue simulator of 
electronic systems, a scenario arising from the intricate relation between statistics and interactions 
in quasi-one-dimensional systems. 
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Experimental developments in ultracold atomic sys- 
tems continue to yield insights into exotic quantum 
phases that were once the dominion of theorists. Since 
the first Bose-Einstein condensates were produced in the 
early nineties P, the steady march of experimental 
progress continues to allow us access to all manner of 
phases from Mott insulator to superfluids [Sj] - see Ref. 
for a recent review. Solid theoretical proposals exist for 
many other unconventional collective states - for exam- 
ple since early proposals to find the elusive supersolid in 
an optical lattice [5| , there has been a spate of recent sug- 
gestions jil-fioj of alternative setups in which a supersolid 
phase may be uncovered. 

A promising possible application of these developments 
is that one might be able to use arrangements of ultracold 
atoms in optical traps as analogue computers for simu- 
lating strong correlations in electronic condensed matter 
systems ll| . As correlations are enhanced by low dimen- 
sions, it becomes natural to consider systems involving 
weakly coupled quasi-one-dimension tubes of particles, 
systems which share their underlying geometry with that 
of some of the more interesting condensed matter sys- 
tems. While it may at first seem preferable to deal with 
fermions if one intends to compare back to electronic sys- 
tems, in this paper we will demonstrate a bosonic model 
with the same low-energy effective theory as a well stud- 
ied model of striped superconductivity in electronic sys- 
tems 12, [13] that has been successfully applied to the 
compound Sr2Cai2Cu2404i [13] but r tiay even have rel- 
evance for high- Tc cuprate materials [l5j. This equiva- 
lence between bosonic and fermionic models arises as a 
result of the quasi-one-dimensional nature of the setup, 
where the role of statistics and interactions are intrin- 
sically intertwined [l6| and implies that suitably chosen 
bosonic experiments can inform us about the low energy 
excitations of fermionic systems. 

Experimental groups now routinely produce such tubes 



FIG. 1: [Color online] Proposed experimental set up. An op- 
tical lattice in the y direction separates the trapped dipolar 
bosons into weakly coupled tubes. The dipole moments of the 
bosons are polarized by an external field, _B, which may be ro- 
tated in three dimensions to change the effective interactions. 

of bosons. The obvious interaction to exploit within 
such systems is the contact interaction as this is tunable 
through Feshbach resonance 



17[. However contact in- 



teractions are too short-ranged to couple distinct tubes 
and so it becomes preferable to deal with the dipole- 
dipole interaction which itself is well tunable in these 
systems [3l and allows the manufacture of custom spa- 
tially anisotropic interactions in such systems |19h22| . 
It is unimportant for the majority of our considerations 
whether the dipole moment is electric or magnetic, how- 
ever one may desire the larger energy scales offered by 
the electric dipole to access some of the physics we will 
encounter. 

In this manuscript we will discuss such a model in 
which a gas of dipolar bosons is arranged in a two- 
dimensional planar array of one-dimensional homoge- 
neous cigars - see Fig.[TJ This may be created by a strong 
two dimensional optical lattice in which the confinement 
in the z-direction completely suppresses hopping in that 
direction, while that in the y-direction allows some weak 
hopping between the tubes. In principle, a weak optical 
lattice may also be present in the longitudinal (x) direc- 
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tion without changing the physics, so long as the lattice 
is incommensurate with the density of bosons. We al- 
low for interactions between the tubes through a dipole- 
dipole interaction which can be tuned via the direction 
of the dipole moment. 

A similar experimental set up has been considered in 
[2^ in the limit of zero inter-tube hopping. The effects of 
this hopping, but in the absence of the dipole-dipole in- 
teraction was discussed in [23|. In this paper, we will con- 
sider both inter-tube couplings to be present, paying par- 
ticular attention to the quantum critical region created 
by their competition. Given the anisotropy of the lattice 
and the freedom of the dipole moment to point anywhere 
in 3-space, we have a great control over the interaction 
strength in the longitudinal and transverse directions. If 
we choose the angle /3 between the dipole-moment and 
the normal to the plane to satisfy sin /? = 1/ {(3 is 
referred to in the literature [21] as the magic angle) then 
interactions in either direction will be repulsive and it is 
this case that we will consider from here on, fixing /3 and 
varying a. 



The dipole-dipole interaction may be split into two 
contributions: intra-tube and inter-tube. While we keep 
explicitly the full structure in the intra-tube part, we ap- 
proximate the inter-tube part by a local nearest-neighbor 
interaction, this approximation is valid for the dipolar in- 
teraction when, as is the case here, the anisotropy of the 
underlying lattice is large. For a fuller discussion see [21j, 
note here though that the dipole interaction in two di- 
mensions (as opposed to three dimensions) is not truly 
long ranged as there is no divergence in the integral of 
the interaction. This means that none of the essential 
physics is changed when considering nearest neighbour 
interactions only as we are regardless operating within 
the correct universality class. Now if x measures along 
tube n (at ?/ = nX) and a is the angle between the hori- 
zontal projection of the dipole moment and the axis of the 
tubes (see Fig. [ij , then the contributions to the dipole- 
dipole interaction take the form 
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dxdx' Pn{x)VlT"{^.^')Pn{x'), 

dxJ± (i't(x)*„+i(x) + *],+i(a;)^'„(a;)) 

dxdx' pn{x)wVlf''{x,x')pn+i{x'). (2) 

Here, ^n{x) denotes the creation field operator in the 
n*^ tube, the hopping parameter J± may be found by in- 
tegrating the confining potential with respect to the ap- 
propriate Wannier functions of the longitudinal lattice, 
and w comes from the integral off the Wannier functions 
in adjacent chains. We may also include a contact in- 
teraction in the above model, however for simplicity we 
consider the case where this is negligible, as in ^255. 

The above described model has two clear strong cou- 
pling limits; when the interaction is dominant, the 
ground state will be a density wave (DW) while if inter- 
actions are weak a superfluid (SF) phase should emerge. 
We will now derive the phase diagram of this model as a 
function of experimentally controllable parameters, pay- 
ing particular attention to the quantum critical region 
between the DW and SF, where we show that a super- 
solid phase may also manifest itself. 

Firstly we bosonize the Hamiltonian according to the 
standard procedure 2^ 2^ via the prescription 
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where po is the particle density (which we assume to be 
the same for each tube), A is a non-universal dimension- 
less constant of order one and ipn and 0„ are the dual 
fields of the tube {[ip{x), dj{x')] = iS{x - x')). 

Now neglecting higher harmonics, additive constants 
and terms which vanish faster than 1/po we arrive at the 
following bosonized description 
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where Cdd is the strength of the dipole-dipole interaction. 



The system is therefore described by a Hamiltonian 
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tions are given by 



^intcr' where the individual contribu- 



-"tube ~ 27r 
P^'^ ~ 27r / a2 



5i„t C0S2((^„ - (pn+l) 



3hop< 



(4) 



Here a ~ pg ^ is the short distance cutoff and the dimen- 
sionless coupling constants gint and ghop are given by 



a'^ PQwCdd cos^ a Ana^AJ^pQ 

5int = , 5hop = — • (5) 
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FIG. 2: [Color online] A contour plot showing the variation of 
the Luttinger parameter K with the dipole angle a and the 
relative strength of the dipole interaction. The most impor- 
tant contour for our purposes is that for A" = ^ as we shall 
see this is where the DW and SF phases are maximally in 
competition. Inset: the bare relation between K and dipole 
strength based on the parameterization in Ref. [2M Although 
for cold atoms with a magnetic dipole moment, poro <gi 1, 
cold molecules with an electric dipole can realistically be in 
the regime where po^'o > 1 29;]. 

In the absence of inter-tube coupling (ghop = 5int — 0) , 
each tube is governed by a Luttinger liquid type Hamil- 
tonian with Luttinger exponent K and sound speed Vg. 
The relationship between these parameters and the orig- 
inal model parameters is non-universal, but has been de- 
termined numerically (28j , where it is found that K varies 
smoothly as a function of po^'o where tq is the effective 
Bohr radius. In our case, tq is given in terms of the dipole 
angle a by ro(a) — ^^sin^a, giving the dependence 
of K on a as shown in Fig. [2] 

This bosonized Hamiltonian is now of a form similar 
to that of a quasi-one-dimensional superconductor, so we 
shall adopt the approach taken in Ref. |l3|. If the in- 
teraction coefficient or the hopping coefficient is clearly 
dominant then in either limit we can consider the mean 
field behavior and arrive at a sine-Gordon model that 
must be solved along with a self-consistency condition. 
Setting /?int — 2i4r^/^ and /3hop — if"^/^, we find a mean 
field Hamiltonian of the form 
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dx 



- 25(cos(/3^)) 



mf 
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Though we explicitly write this for the interaction 
dominated case, due to the duality between 9 and ip we 
can treat both regimes simultaneously. 

The expectation value (cos(/3(^)) is known within the 
sine-Gordon model [s^ so we can solve the self consis- 
tency conditions to find that g™^ ~ 2(72(1+1^)^ where 



d = 13"^/ 4: is the scaling dimension of the operator in 
question. The scaling dimensions of our two inter-chain 
perturbations are di^t = K and c?hop = ^/{'^K). 

Inspecting the mean field solutions we have acquired, 
we see that the hopping dominated region is character- 
ized by a non-zero expectation of {cos{9)) which cor- 
responds to ^ 0, and as such represents a 
SF phase. In the opposite case when the interaction 
dominates, the phase is characterized by a non-zero ex- 
pectation of (cos(2i^)) which corresponds to {pn{x)) = 
Po + C(— 1)" cos(27rpoa;), and hence represents a checker- 
board DW phase. Notice that if we were to abandon the 
"magic angle" which is holding all interactions repulsive 
and allow for attractive interactions, the DW phase can 
also be driven into a striped phase, a detailed study of the 
phase boundary of these two DW phases is undertaken 
in [2^. By changing /3 one could see either DW phase 
competing with the SF phase, for concreteness we will 
hereon continue to assume that the relevant DW phase 
is the checkerboard. 

In order to determine which phase is dominant for a 
given parameter range, and hence find the phase diagram, 
we consider both the hopping and interaction terms as 
perturbations over a Luttinger liquid and ask which has 
the higher critical temperature. We find the SF and DW 
susceptibilities (xsf and Xdw respectively) within the ran- 
dom phase approximation using x — X^*'V(1 ~ ijrxf*'^)- 
The LL susceptibility for either phase is (see e.g. [27|) 
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The critical temperature Tc is estimated as the tem- 
perature at which an instability arises from a diver- 
gence in the associated susceptibility. As our system is 
two-dimensional, these transitions are in the Berezinskii- 
Kosterlitz-Thouless universality class, however this sim- 
ple calculation still estimates correctly 13| in the weak 
coupling regime. Hence we find Tc from the condition 
-^/(o) ^ -2_ so we find that the interactions dominate and 

1 2K 

the system takes on a DW phase when 3;^^^'^ ^ 9hop^ ■ 
In the opposing region the inter-site hopping dominates 
and the system becomes a superfluid. 

The parameter we wish to focus upon is the angle a 
which comes from the geometry of our proposed setup. 
The explicit dependence of the critical temperature of 
the SF and DW phases on this geometric parameter are 
given by 



ksT^^ — ^ (^Shop 



iK[Rgf^^f^ Bin- 



iff-rcos^a) 



(8) 
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FIG. 3: [Color online] A plot of the critical temperature for 
the DW and SF phases and the implied phase diagram for 
an illustrative set of parameters (gint = 0.1, g^op = 0.2 and 
R — 250) chosen such that K ^ ^ at the critical point. Inset: 
correction the phase diagram at the critical angle showing 
the onset of a supersolid (SS) phase. Based on realistic ex- 
perimental parameters for the dipole moment of e.g. KRb 
[33 | and a lattice spacing A ~ 500nm, we estimate the max- 
imum critical temperature of the DW phase to be of order 
Tdw ~ 50nK. 



where R — 



andgfnf'' = niax(gint(a)) = gi„t(a = 



0). Fig. [3] shows the dependence of these critical temper- 
atures on a for indicative values of R, (/hop and g™^^. 

The above quasi-one-dimensional analysis is valid as 
long as the parameters are sufficiently far from critical- 
ity. In fact, even though the ordered phases are necessar- 
ily two-dimensional, deep in either of them, the correla- 
tion functions retain their one-dimensional character 31 1 . 



However, near the quantum-critical point (QCP) between 
the DW and SF, the entire quasi-one-dimensional formal- 
ism breaks down, and a full two-dimensional field theory 
is needed, the dimensional anisotropy here being an ir- 
relevant operator [34 1. 

Exactly at the QCP, where DW can be smoothly 
rotated into SF, the order parameter acquires an en- 
larged (and non-Abelian) symmetry [l3l |. This is in 
analogy with SO{5) theories of the antiferromagnetic- 
superconductor transition [s^], except that the enhanced 
symmetry at the SF-DW transition is 0(4)[3^. Follow- 
ing Jaefari et al. [s^], in the vicinity of the QCP we can 
then write a 2 -I- 1 dimensional non-linear sigma model in 
terms of this enhanced 0(4) order parameter, which has 
a Lagrangian 



1 



^50 



wd^NaOabdf^m - hNaOabNb. 

(9) 

Here, TV is the order parameter, go ^ \/gint + ffhop is the 
effective coupling constant, w ^ K — signifies the de- 
parture from the point where the two inter-chain terms 
have the same scaling dimension, h gi^t — 5hop mea- 



sures the distance from the critical point, and O and O 
are two fixed matrices that break the enhanced symmetry 
in the appropriate way. It turns out [s^l that the term 
proportional to w is marginally irrelevant in the renor- 
malization group (RG) sense, meaning that the generic 
quantum critical behavior of our original model may be 
described by the above non-linear sigma model. Note 
that all we are doing here is writing the most general 
Langrangian with all the correct symmetries. This model 
has been well studied using RG methods 34 , 3^ where it 
is found that there is a tetracritical point at (70 = 5c 1 
separating different types of critical behavior. In the case 
go < 5c which is always satisfied in our present weakly- 
coupled chains scenario, there is a region around ttc where 
there is coexistence of both the DW and SF order param- 
eters, in other words a super-solid phase. The corrected 
phase diagram in the vicinity of the critical point is shown 
schematically in the inset to Fig. [S] 

Unfortunately precise predictions regarding the super- 
solid phase are not easily attainable because this phase is 
deep within the strongly correlated regime and not acces- 
sible from our interchain perturbation theories. Its exis- 
tence must rather be inferred based on RG relating to the 
symmetries of the order parameter at the quantum crit- 
ical point, leading to a qualitative phase diagram. This 
said, as indicated in to Fig. [3J the relevant region to look 
for the supersolid phase in any experimental investiga- 
tions of this system is set by the transition temperature 
of the suppressed phase, which could be found by our 
mean field calculations. 

We can compare our results to those found in isotropic 
2D analogs of this system: in an isotropic lattice with 
nearest neighbour interactions it has been shown that Q 
the supersolid phase may arise at the interface of the SF 
and DW phases. The stability of the SS phase is consid- 
ered in [36| in which they note that this phase is possible 
in the presence of a lattice. This is as opposed to the free 
system in which, for a dipole-dipole interaction, it has 
been seen [s^l that there is a depletion of the SF phase 
before the onset of a dipole solid phase. Although our 
method is not applicable to the isotropic limit, the fact 
that the phase diagrams for models on a 2D lattice are 
so similar suggests that there will not be any additional 
transitions as we decrease the anisotropy. 

Returning finally to the proposed experimental setup 
in which we expect this phase diagram to manifest itself; 
some groups [38| have reported single-site resolution in 
fiourescence imaging of ^^Rb in a two dimensional lat- 
tice setup similar to that we describe so it may well be 
possible to probe the DW phase of this system in situ. 
A common way to probe the SF phase is releasing the 
bosons from their confinement, and looking for high con- 
trast in the interference pattern in the absorbtion data 
[sj. The hallmark of the supersolid phase we have de- 
scribed would be the presence of both a DW and SF 
order parameter. 
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In summary, we have shown how a popular model in 
condensed matter physics of striped superconductivity 
may also appear in a cold atom experiment with trapped 
dipolar bosons in a quasi-one-dimensional geometry. The 
rich phase diagram of the model may be easily explored 
by changing the angle of a polarizing field, and shows 
crystalline, superfluid and supersolid phases, on top of 
the high-temperature sliding Luttinger liquid phase. 

We thank Joseph Betouras, Dima Gangardt, Chris 
Hooley and Jorge Quintanilla for useful discussions. 
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